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1 Introduction 



A very convenient way to non-perturbatively regularize a QFT is to put the dynamical vari- 
ables of the theory on a regular spacetime lattice (in the functional-integral formulation) or 
on a regular space lattice (in the hamiltonian framework). This introduces a "natural" cutoff, 
roughtly equal to the inverse of the lattice spacing, either on both energy and momentum or on 
momentum alone. Usually this procedure breaks the symmetry properties of the action down to 
a lower level: Lorentz or Euclidean invariance reduces to invariance under discrete subgroups, 
scale invariance in massless theory is broken explicitly by the cutoff, and very often also internal 
symmetries, either global or local, are difficult to keep. 

Therefore it is very interesting to find regularization procedures that preserve as much as 
possible of the characteristics of the continuum theory. This issue is particularly important in 
the case of two-dimensional models which are integrable at tree level and are supposed to be 
so also at the full quantum level. One would like to have a non-perturbative lattice definition 
of such quantum theories which preserves integrability. 

A quite general solution to this problem is based on the so-called light-cone approach || , in 
which the 2D Minkowski spacetime is discretized in light-cone coordinates. The basic object in 
this approach is the 72-matrix, that is a solution of the Yang-Baxter equations which character- 
ize the factorized scattering of a 2D integrable QFT. This i?-matrix is regarded as a collection 
of quantum amplitudes for the scattering of "bare" objects, which move with the rapidity cutoff 
0, on each vertex of the light-cone lattice, casting the model in question in the form of a vertex 
model. Then the full machinery based on monodromy and transfer matrices WO Q can be 
set up and the algebrized or analityc Bethe ansatz (BA) can be used to completely diagonalize 
the transfer matrix and, with it, the total momentum, the Hamiltonian and all other conserved 
charges. The continuum limit may then be explicitly performed by letting O go to infinity in a 
well defined way as the lattice spacing vanishes. 

A drawback of the standard light-cone approach is the nonlocality of the lattice Hamiltonian. 
While this does not constitute a real problem for the continuum limit, either at the bare or 
renormalized level, it makes more difficult to properly handle the full excitation spectrum 
and to study the conformal limit, which allows to identify the integrable model at hand as a 
perturbed CFT. A sligthly modified version of the light-cone approach without such difficulties 
was recently put forward in Q: rather than as logarithm of the unit time evolution operator (or 
diagonal-to-diagonal transfer matrix), the lattice Hamiltonian is identified as the first of the 
series of local charges obtained by suitably differentiating the alternating transfer matrix with 
respect to the spectral parameter. Such identification was made before, whithin a different 
context, in The basic property of this modified approach is the locality of the lattice 

Hamiltonian, which allows to safely regard the time as continuous while the space is still discrete, 
restricting the UV cutoff only to the space momentum. 

In this work we present a detailed application of the local light-cone approach to the massive 
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Thirring model. This is probably the simplest case, being based on the well known, almost 
paradigmatic six-vertex i?-matrix, without any quantum group restriction, and was for such 
reason the first model studied also in the nonlocal approach Q. Nonetheless there are some 
interesting non-trivial points that require a careful examination. 

First of all one must take into account the Nielsen-Ninomiya theorem [[H]] , since the lattice 
Hamiltonian is local and chiral-invariant in the — > oo limit (this is one of the most important 
differences of the light-cone approach with respect to Liischer's regularization based on the 
XYZ spin chain |l3|]: the latter is indeed integrable but has neither U(l) invariance nor a local 
implementation of chiral transformations). As a consequence one finds the "fermion doublers" 
both in the perturbative spectrum and in the exact Bethe ansatz spectrum. It is then important 
to check whether these massless doublers indeed decouple from the massive Thirring particles. 
We show the answer to be affirmative even off shell, for the local continuum fields, although 
the mechanism is quite non-trivial. 

Secondarily, we examine in detail the problem of boundary conditions and their effects on 
the exact spectrum. In particular, by carefully handling a completely fermionic formulation we 
are able to show that the excitations over the ground state carry the correct U(l) charge which 
corresponds to dressed fermions interpolated by the bare fields. This should be compared with 
the result proper of the periodic spin chain, with excitations carrying half the U(l) charge of 
the fermions. 

Another interesting point concerns the structure of the perturbative vacuum on the lattice: 
while the one-particle spectrum over the emptied Dirac sea (the state killed by the local fermi 
fields) has a anisotropy-dependent zeroes and no simmetry between positive and negative ener- 
gies, this simmetry is restored and the anomalous zeroes move to the boundary of the Brillouin 
zone simply by normal-ordering the U(l) currents in the lattice Hamiltonian. This facts allows 
to isolate the effects of the interaction, even before the continuum limit, in a cutoff-dependent 
mass renormalization and in a finite rescaling of the velocity of light. 

The finite renormalization of the speed of light is one last subtlety that requires a proper 
treatment. While such renormalization is absent in the nonlocal light -cone approach, where 
the simmetry between space and time is mantained all along, nothing forbids it in the local 
formulation, since time may be regarded as already continuous while space is still discrete. We 
handle this by intruducing a time unit at which is independent from the lattice spacing a of 
the space chain. The velocity of light, either bare or renormalized, emerges quite naturally as 
finite ratio a/at- 

This paper is organized as follows. In section || we describe the basic framework of vertex 
models and derive in a purely algebraic way the local lattice Hamiltonian, using first the R— 
matrix written in spin language. In section |3] we discuss the subleties related to the formulation 
on the ligh-cone lattice of the system using a fermionic approach. Indeed the translation of the 
i?-matrix in fermionic variable is quick (after taking in account some important changes of sign 
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due to the Fermi statistic), but involves a careful definition of the boundary condition. The 
explicit form of the hamiltonian and the dispersion laws for the lattice fermions derived in this 
fermionic setup are shown in section ||, where it is also discussed the normal-ordering prescrip- 
tion we adopt for the U(l) currents over a completely occupied Dirac sea. This corresponds 
to an antiferromagnetic ground state in spin language. The continuum limit is considered in 
section |5[ where abelian bosonization tricks are used to disentangle the mixed currents terms 
that arise in the naive continuum limit. In this way we shows that in the continuum Hamilto- 
nian does describe two fermi fields, one massless and one massive. In section ^ the results of 
the Bethe ansatz are briefly rewieved, showing some novelty regarding the meaning of the hole 
charge in the framework with antiperiodic boundary conditions and the matching between the 
dispersion laws perturbatively derived from the lattice Hamiltonian and the exact one based on 
the Bethe ansatz. Finally, in [?], we study the effects of renormalization and of the trasforma- 
tion to the decoupled description on the relation between the various coupling constants: for 
instance, the current-current Thirring coupling and the sine-Gordon coupling constant f3 are 
related in the standard one only after a suitable power serie redefinition. Some comments on 
the results obtained and on possible further developments can be found in M. 



2 The basic framework 

It is well known ||] || that the 6V model, as well as the XXZ spin chain related to it, may 
be formulated starting from a collection of two-dimensional vector spaces {Vj, j = 1,2 ... , N} 
and local R-matrices Rij acting on the tensorial product Vi (8> Vj of two such spaces. These Rij 
are written in terms of the Pauli's matrices crj, crj, cr|, j = 1, 2, . . . , N, as 

*y (A) = ^ + HT^J + b [°t°i + °T°j] 

where the = a x ± ia y and the trigonometric Boltzmann weights b, c are parametrized as 
follows by the spectral parameter A: 

sinh A 



6(A) 



sinh (i'y — A) 



i sm 7 

C = c(A)5 sinh(H-A) < 2 - 2) 

The choice of weights made here guarantees that the .R-matrices are unitary for real A and 7, 
that is B}- k Rjk = 1. It is straightforward to check that this reduces to the identities |6| 2 + |c| 2 = 1 
and bc + bc = 0. As we shall see below, the unitarity property is important in order to interpret 
the transfer matrix as a temporal evolution operator. The regularity condition of the R-matrix 
is fulfilled by eq.(|2.1j) as Rjk(0) = 1. Most importantly, the i?-matrices satisfy the Yang-Baxter 
equations (YBE)0] 

R ij (X)R jk (\ + ^Rijin) = R jk {fi)Rij{\ + fi)Rjk(X) (2.3) 
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which ensure the integrability of the 6V model in any framework. 
The 'bare scattering' 5-matrices are defined as 

Sij(X) = P ij R ij (X) (2.4) 

where the permutation operators Pij interchange the vector space Vi and Vj. PijVi®Vj = Vj<8>Vj. 
In terms of the 5-matrices the YBE eq.(|2.3|) can be reformulated as 



Sjk(X)S ik (X + n)Sij(fi) = Sij(fi)S ik (X + n)Sjk(X) (2.5) 

Let's now introduce the fully inhomogeneous monodromy matrix T(A|{#j}) associated with 
the auxiliary "horizontal" vector space Vo 

T(X\{9 i }) = S 10 (X + 9 1 )S 20 (X + 9 2 )...S N o(X + e N )=(t (2.6) 

where the operators A, B, C, D act in the full Hilbert space Vi <8> V2 • ■ • (8> Vn- The monodromy 
matrix, thanks to the YBE, satisfies the Yang-Baxter algebra (YBA) 

R(X - [T(X\ {9 { }) ® {e i })} = [T(v\ ® T(X\ {6 t })\ R(X - /i) . (2.7) 

These implies a set of commutation rules for yl, C, D, among which the following play a 
central role in the algebraic Bethe ansatz: 

b(fi- X)A(X)B(fi) = +B(fi)A(X)-c(fi-X)B(X)A(fi) 
b(\-n)D(\)B(ji) = +B(fi)D(X) - c(X- fi)B(X)D(fi) (2.8) 
B(X)B(p) = B(p)B(X) . 

Taking the trace of the monodromy matrix over the horizontal space we obtain the transfer 
matrix 

t(X\{9 i })=tv T(X\{9 i }) (2.9) 

For fixed arbitrary set of vertical inhomogeneities thanks again to the YBE, the transfer 
matrices form an infinite set of commuting operators. 

[t(x\{e i }),t{ Ji \{e i })] = o 

Since we are trying to regolarize a relativistic QFT on a light-cone lattice, we choose the vertical 
inhomogeneities in a particular way, consistent with the propagation of 'bare particles' moving 
along the two diagonal directions with cutoff rapidity ±0, respectively: 

^ = (-l) m 9, i = l,2,...,2N . (2.10) 

We have changed N to 2N to ensure periodic boundary conditions. 
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Inserting these alternating inhomogeneities in eq.([2.6|) we obtain the alternating monodromy 
matrix 

T(X\ 9) = 5 M (A + O)S 20 (X - 6) . . . S 2N0 (X - 6) (2.11) 

and, taking the trace as in eq.( |2.9| , the alternating transfer matrix t(X\ 0) = troT(A| 0). 
The regularity condition Rjk(0) = 1 and the permutation algebra 

{A-knPij h J) k, n all distinct 
AnPij j = k; i, j, n all distinct , (2.12) 
j = re; i, j, k all distinct 

which holds for any operator Aij acting nontrivially only on Vi (8> Vj, imply the fundamental 
relation 

t(Q\ Q) = U L , t{-Q\ 9) = 4 • (2-13) 

Here Ur and Ul are the right and left diagonal transfer matrices (they move by one lattice 
spacing in right-upward x + t and left-upward x — t direction respectively) defined as 

U L = VR 12 Ru . . . R 2N -i 2N (2.14) 
U R = V^RuRu . . . R 2N -x 2N (2.15) 

where Rjk = Rjk(2@) and V is the left shift operator V = P\ 2 nP22N ■ ■ ■ -fW-i27V- The 
derivation of these formulae is purely algebraic; for Ul(@) we have 

N 

t(Q\Q) = troll S 2 j-io(2e)P 2j0 

j=l 

~N-1 

= (tr P 2 N ) J\ S 2 j-i 2 N(2Q)P 2 j 2 N 5 , 2AT_l2Ar(20) 

J= l 

f N-l \ N 

II P ^2n \ Y[P 2j -i 2j R 2j -i 2j (2Q) 

J=l } 3=1 

N 

= V]J i? 2i _i 2i (2G) 

3=1 

= U L (2.16) 

with a similar calculation for Ur. 

The unit time evolution operator is U = UrIIl- it causes a displacement at, the lattice spac- 
ing in the time direction, upwards on the light-cone lattice, leading to the following definition 
of the lattice Hamiltonian: 

H = iat 1 logU. (2.17) 

Evidently this Hamiltonian is nonlocal. Similarly nonlocal is the lattice momentum operator, 
naturally defined as 

P = -ia' 1 logV 2 , 
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where a is the lattice spacing in the space direction. On the other hand, from the commuting 
family of alternating transfer matrices it is possible to obtain a full hierarchy of local charges 
in involution. It suffices to take the logaritmic derivative of t(A| O) with respect to the spectral 
parameter at A = ±0: 



H, 



± 



-l-ra 



Qn 



logi(A|G) 



A=±e 



(2.18) 



By purely algebraic calculations similar to those of eq.( [2.16D , one verifies that H^(Q) couples 
2n + 1 neighboring sites. Unlike in eq.( |2.16| ), in this derivation it is crucial that the i?-matrices 
satisfy the YBE. Since H^(Q) commutes also with t/(0), it is a conserved charge. 
The charges of level 1 read 



N N 

H i = E /j 2 i -i (20) , #r = E h v 1 

3=1 3=1 



-20) 



in terms of the 'Hamiltonian density' 

KiW = — -Rnn+l(A) 1 R nn +l(\) + R n -i n (0)R nn+ i(\) 

With them, one can now define the local Hamiltonian 



H = 

2a t 



H+ + Hi 



(2.19) 



(2.20) 



(2.21) 



which is indeed hermitean thanks to the unitarity of R— matrix. Of course, with this choice of 
hamiltonian, the evolution operator is U(t) = e~ ltH , with the time t continuous and at merely 
fixing the scale of time or energy. 



3 Fermionic formulation 

The U(l) invariance of the 6V R— matrix corresponds, in the light-cone framework, to the 
conservation of bare particles. In fact the ferromagnetic state with all spins up, | + + . . . +), 
may be regarded as 'bare vacuum state' (the state with no bare particles). Then we can say 
that a state with r flipped spins located at 1 < ji < j'2 . . . < j r < 2./V, that is the state 

\ji,h, ■ ■ ■ Jr) = vj^Jz ■ ■ -<rj r |+ + • • •+) 

contains exactly r bare particles at the same locations. The particle number r is conserved in 
time, that is along the vertical direction throughout the lattice, thanks to the U(l) invariance 
of the R— matrix. 

It is clear that these particles are identical and satisfy the Pauli exclusion principle, since 
(crj") 2 = 0. On the other hand, since \aj , a^] = for j ^ k, they are of bosonic type. This can 
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be remedied by means of the well known Jordan- Wigner transformation from the spin operators 



Oj and (Tj to lattice fermion fields 



j-l j'-i 

^•=<n< ^j=^n< (3.1) 



n=l n=l 



satisfying the canonical anticommutation rules 

, a} = o = {vj , 4} Hi > 4> = ^ • ( 3 - 2 ) 



The string of in eq.( ]3.1| ) has a nontrivial effect only on the boundary conditions. In fact 
it cancels completely out of all local R— matrices with neighboring indices such Rjj+i (with 
1 < j < 2N — 1), which have the fermionic form 

R jj+1 = 1 + bK jj+x + (c - l)(Qj - Q J+ if (3.3) 

where 



Kij = + 1>tyi , Qj = 4>fy = |(i - of) • 

Thus the string of <r^ would also drop out of the evolution operator U and of the local hamilto- 
nian Hi, if it were not for the periodic boundary conditions. The troblesome object is R2N i(20), 
which reads in terms of fermion operators 

R 2N1 = 1 - b [4 N M-) F + (-) F 4^2n] + (c - l)(Qi - Q2N? (3.4) 

where 

27V 

(-f = n °* = e,7rQ 
3=1 

is the longest possible string, that is the fermion signature, and Q = J2j Qj is the total bare 
particle number. Similarly, since the left shift operator V acts on the Pauli matrices as 

V^CTjV = <Xj + l 

it cannot shift exactly also the fermion fields. Rather we have 

VHiV = fl < = ^ + ie^ , (3.5) 

n=2 

for j = 1,2,... ,2N- 1, and 

V^2NV = 1p2N+l = -lpl(-) F ■ 

Together with eq.(3.4), this last relation suggests that PBC on the spin operators become a 
sort of F-twisted boundary conditions 

i>2N+i = -M-f (3-6) 
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on the fermion operators. However, if it is true that this guarantees R2N2N+1 = R2N1, eq.(jO|) 
prevents the identification of V with an exponential of the fermion total momentum. In par- 
ticular, V 2N is the identity in the full vector space V (27V) , and hence V- 2N ipjV 2N = xpj, which 
shows the conflit between ipj+2N = V~ 2N ipjV 2N and the extension tpj+2N = ~i ) j(~) F °f the 
F-twisted relation ( |3.6| ) to all fermion operators. Of course, we could define the true fermion 
shift operator V through 

for any j, and impose uniform F-twisted boundary conditions via 

V~ 2N ^V 2N = -^j(-) F . (3.7) 

Certainly V commutes with U and Hi but, unlike them, it is not related in any obvious way to 
the transfer matrix t(X\ 0), which is the object we are able to actually diagonalize by means of 
the algebraic BA. Therefore the translation of the light -cone 6V model and its BA solution from 
its original spin formulation into a fermionic theory, by means of a straightforward application 
of the Jordan-Wigner transformation, remains unsatisfactory due to boundary effects. 

Although we expect that these boundary effects will loose importance in the limit N — > 00, 
it is convenient to look for a purely fermionic formulation, in which all basic objects, like 
R— matrices and exchange operators, are written from the start in term of fermion fields for 
any pair of indices. 

To this end, let us notice that the matrices Rjj+\(X), whether written in spin (eq.( |2.1[) ) or 
fermion language (eq.( |3.3| )), satisfy the YBE in the restricted form 

j+ift + VjRj-ijW = Rj j+i{^)Rj-\j{^ + ^)Rj j+iW ■ 

But since j — 1, j and j + 1 simply refer to three distinct anticommuting fermions, the matrices 

Rij(X) = 1 + b(X)K i:j + [c(A) - l](Qi -Qj) 2 

will fulfill the general form ( [2. 3D of the YBE, providing another solution distinct from Rij(X). 
In fact, Rij(X) / Rij(X) for \i — j\ > 1. 

Next we build the fermion permutation operators Pij, defined by the relations 

They are written in terms of the fields simply as 

Pij — 1 -\- Kij . 

Then we can build the S— matrices 

S tj = PyRij = 1 - 2Q i Q j + cKij + (b- l)(Qi - Qj) 2 . (3.8) 
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Unlike in the spin framework, now the relation between 6V R— and 5— matrix does not reduce 
simply to the exchange c # b, since Fermi statistics requires that must be —1 in the doubly 
occupied state, rather than 1. This is taken care by the last term in eq.(3.8). The matrices Sij 
and Rij manifestly commute with the bare particle number Q which generates the symmetry 
group U(l). To lightens the notation, from now on we drop the "throughout, reinstating it only 
when strictly necessary. 

We have now all the ingredients to build the relevant global objects, which are the fermionic 
analog of V, Ur, Ul, T(X\ {9%}) and t(A| {Oi}), with all the relations that we found in section (|2| 
valid also for the new objects, since they are based solely on algebraic properties like regularity, 
YB algebra and permutation algebra. In particular, the alternating monodromy matrix 

T = T(X\ {9i}) = S10S20 ■ ■ ■ S2N0 

can be written 

T = A + B^ Q + C^l + (D- A)ip% 

where ipo and ip\ are new auxiliary fermion operators anticommuting with all the previous ones, 
and A, B, C, D are global operators in the full fermionic Fock space. Notice that ipQ commutes 
with A and D but anticommutes with B and C. In fact, one easily verifies that A and D have 
an even fermionic grade (that is they are sums of terms containing an even number of ipj and 
ipj, j = 1, . . . , 2N), while B and C have an odd fermionic grade. 

To write the YB algebra it is convenient to rename ipo into, say, xi-, an d introduce another 
pair X2, X%i anticommuting with all ipj, j = 1, . . . , 2N as well as with xi- Then we can write 

T r = A + B Xr + Cxt + iP- A)xtxr 

and the YB algebra takes the form of Eq. ( |2.7| ) 

Su(X - fi) Ti(A| {Oi}) T 2 (jjl\ {e % }) = r 2 ( M | {Oi}) Ti(A| {Oi}) S 12 (X - fj) 

where (see eq.plf)) 

S12 = 1 + c[x\x2 + X2X1) + {b- + X2X2) ~ 26x1X1X2X2 

To obtain the commutation rules for A, B, C, D the algebra is now straightforward: one finds 
some differences of sign with respect to the rules expressed in eq.(|2.8|), namely 



b(jji-X)A(X)B(jM) = +B((j,)A(X) - c(n- X)B{X)A(fi) 
b(X- fi)D(X)B{fi) = -B(n)D(X) + c(X- fi)B(X)D{fi) (3.9) 
B(X)B(jm) = -B(fi)B(X) . 

Of course, the anticommuting nature of the "creation operators" B(X) appears very natural in 
this fermionic setup. The other changes of sign concern only the commutation rules between 
-D(A) and B({i), and could be traced to the fact that Sij = — 1 in the doubly occupied state. 
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One last subtlety concerns the meaning of the trace operation. In this fermionic setup the 
correct definition would be 

t = tr T = (0|T|0> — (1| T |1) = A — D 

where |0) is the state with no auxiliary fermion and |1) the state with one. Then tro Pjo = 1 
for any j and we find the fermionic light-cone version of eqs. (|2.13| ) in the form 

t(&\ Q) = U L , t(-@\ 0) = £/£ . (3.10) 

This choice corresponds to periodic boundary conditions on the fermions, that is ^j+2N = ipj- 
On the other hand, we may take as trace what is other contexts is actually called 'supertrace', 
that is 

t' = str T = (0|T|0) + (1|T|1) = A + D . 
Then we find stro Pjo = 1 — 2Qj and correspondingly (see eqs.( |2.16| and eq.( |2.13| )) 

t'(@\ 0) = (1 - 2Q 2N )U L , t'(-e\ 0) = 4(1 - 2Qi) . 
In this case we could take the unit-time evolution operator to be 

e -iaH = e)t'(-e| ©) f = u' L u R 



where as before U\ = R12R34 ■ ■ ■ R2N-12N (see eq. (|2.14|) ), while 

U' L = (l-2Q 2N )U L (l-2Q 1 ) 
= (l-2Qi)C/ 2 (l-2Q 2A r) 

= ^23^45---(l-2Qi)i?27Vi(l-2Qi)y . (3.11) 

Similarly we now define the unit-space traslation as 

e iaP = t'(e\ 0)t'(-0| 0) = [(1 - 2Q 2N )V] 2 . 

Since (1 — 2Qj)ipj(l — 2Qj) = ipj, this choices correspond to antiperiodic b.c. on the fermion 
fields 

1P3+2N = e- lLP ^e iLP = (3.12) 

where we have introduced the spatial size of the system L = Na. 

In summary, we see that for both choices of trace, leading to either periodic or antiperiodic 
fermions, as well as in the case of periodic spins, the nonlocal hamiltonian and total momentum 
are related to the tranfer matrix as 

e-™ 6 =i(0|O)t(-0|O)t , e iaP = i(0| O)t(-0| 0) (3.13) 

where we may now drop the ' for the antiperiodic case, provided we keep in mind the two 
different ways in which t(X\ 0) is written in terms of the diagonal elements of the monodromy 
matrix, either A — D or A + D. It should be clear that identical conclusions about the the b.c. 
apply in the framework based on the local hamiltonian of Eq. ( [2.21D . 



10 



4 Explicit form of the local hamiltonian 



We shall now obtain the explicit form of the local hamiltonian H in terms of the fermionic fields 
ipj, j = 1,2, ... , 2N. By means of the Jordan- Wigner transformation one can always revert to 
the spin formulation, keeping in mind the effects on the boundary conditions. For definiteness 
we shall choose the antiperiodic b.c. for the fermion fields. When we insert the expression ( p.8| ) 
for the 6V R— matrix into the formula for the hamiltonian density 6(A) (see eq.(2.19)), we need 
to perform rather long albeit trivial algebraic manipulations with the fermi fields. In particular 
we find 



Rj n (\y Rj n (X) 

Rj n (X) Jf Rij(0)Rj n (X) 



= (be + bc)K jn + (66 + cc)K 2 jn 

= 6 [4(^n + c^)+h.c] 

+ 6o[bcipjip n + h.c.) + Qi + ccQj + bbQ n ] 

+ 6 [(6 + b)QjK in + (c - c)Q n (iptyi - ^Vj)] 

- 2c [bcQi(i(j\ij n - iplipj) + Qi(ccQj + bbQn)} 



(4.1) 



(4.2) 



Tit rriTjl ~ ntKy^ysj 

where b = 6(A), c = c(A), 6o = 6'(0) and Co = c'(0). In the derivation of these results the 
unitarity relations bb + cc\ = 1 and be + be = have been used. To obtain H we must now set 
(i,j,n) = (j — l,j, j + 1), then put A = 26 when j is odd and A = —20 when j is even, and 
finally sum up over j. H is the sum of a piece quadratic in the fields and a piece quartic in 
them 



H 

H 2 



«2 + H4 

-I N 



(4.3) 



-a* 



2 sin 7 



i=i 



H 4 



+ 2(v + cos 7 )(Q 2i _i + Q2j)] 
„-i N 

{~(b + b){Q2j-\K 2 j-22j + Q2jK 2 j-12j+l) 



2 sin 7 



where 



+ (C - C) [Q2j(lp 2j ^2j-2 ~ Vlj-2^2i-l) + Q2j+l(lp 2j ^2j ~ ^ 2j ^2j-l) 
+ 2lW COS 7 Q 2 j-2W 2 j-X^2j ~ Ipljfoj-l) + Q2j-l(lp 2j+ i^2j ~ Ipljfoj+l) 

- 2(v + cos7«)Q 2 j (Q2j-i + <?2j+l) - 2cos7&6((5 2 j-iQ2i+i + Q2.7Q2.J+2)} 

it = i sin 7(60 + &c) 

v = i(bb + cc) 

w = ibc = —ibc 

d = u + iw cos 7 + c , 



11 



The quadratic part H2 is better analyzed via the following Fourier transformation (recall the 
antiper iodic b.c.) 



1 



- y 



-n<q< 



2vr 
9G iV 



(4.4) 



In the limit N — > 00 of an infinite chain, the sum over g becomes an integral over q running in 
the first Brillouin zone (— vr,7r). Then H2 takes the form 



H 2 = r dq^h{q)^{q) 

J —IT 



where h(q) is the two-by-two matrix 



h(q) 



a^ 1 ( be~ iq + be iq + 2(v + cos 7) 



2 sin 7 



(2(1 + e^) 



d(l + e~ iq ) 
be iq + be~ iq + 2(v + C0S7) 



The two eigenvalues of h(q) represent the bare energy branches of our lattice model 

E±{q) = [2{v + cos 7) + (b + b) cos q±([-(b- If sin 2 g + 2dd(l + cosg)] 1 / 2 } (4.5) 

where C is +1 in the first Brillouin zone and in all the odd ones, while it is —1 in the even zones. 
These dispersion relations are depicted in fig.l. 



-0.5 




-8 -6 -4 -2 2 4 6 8 



Fig.l: Energy branches for 6 = 2 and 7 = Qn/W in units of a. 



Evidently all negative energy levels, for both branches within the first Brillouin zone, should 
be filled to obtain the lowest energy state. One must take into account, however, that these 
bare fermions are interacting and that this might very well change the shape of the dispersion 
relations themselves. The algebraic Bethe Ansatz will take care of this exactly. At this stage it 
is enough to assume, as natural, that in the interaction picture there exist an equal amount of 
positive and negative energy levels, so that the perturbative filled Dirac sea (the perturbative 
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vacuum state of the QFT) is characterized by half-filling, namely (Qj) = 1/2. This is an 
antiferromagnetic state in spin language. It giustifies the following normal-ordering prescription 

1 



Qn — • Qn ■ ~t~~ 



(4.6) 



which has a dramatic effect on the quadratic part of the Hamiltonian, leading to the perturbative 
one-particle energy spectrum (see fig.2) 

1/2 



E = ±\Qai l 



sinh 40 



sinh 20 + sin 7 



sin 2 q + -(moaj) 2 (l + cos q) 



where 



m = 2a' t 



_ 1 sin 7 e^oo 



sinh(20) 



4a t 1 sin 7 e 26 



(4.7) 
(4.8) 




Fig.2: Energy branches after normal-ordering for 9 = 2 and 7 = 67r/10 in units of a t 



This dispersion relations are manifestly simmetric under reversal of energy, showing the self- 
consistency of our normal-ordering assumption. Once all negative energies in (fO|) are filled, 
one obtains a positive spectrum of particles and holes all with the positive energy of eq.( |4.7l ). 

It is also clear that eq.(|4.7|) represents a lattice approximation to the relativistic spectrum 
of massive particles. To see this we set q = pa and let a, at — * 0. Then we obtain 



E 



cq\Jp 2 + mfcl 



where Co = a/a t is the (bare) velocity of light. It appears natural to choose spacetime units so 
that Co = 1 . Of course one should expect this choice not to necessarily work in the renormalized 
limit to be discussed later. 



The dispersion laws of eq.(4.7) has a peculiarity though: it also describes massless particles 
at the boundaries of the first Brillouin zone. This is inevitable, since we are working with a local 
lattice Hamiltonian which for O — * 00, that is in the massless limit too — > 0, becomes chiral 
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invariant. The Nielsen-Ninomiya theorem [jTo| then implies the existence of the (in) famous 
'fermion doublers'. In the model at hand, the left and right modes around q = are massive for 
finite 0, while the left and right doubler around q = ±7r remain massless. In the limit — > oo 
at fixed lattice spacing the model becomes gapless and it corresponds therefore to a regularized 
Conformal Field Theory. According to the general rules, the neighborood of the critical point 
= oo defines a regularized Perturbed CFT. By letting — > oo and a — > simultanously in a 
suitable way one recovers the continuum PCFT. The CFT describing the critical point and the 
perturbing operator will be identified in the next section. 



5 The continuum limit 

We now consider the continuum limit a — ► where only the small energy excitations (as 
compared with a™ 1 ) of the fields are retained and the massless dispersion relations are linearized 
around their zeroes Q . In this limit the bare mass mo is kept fixed (the renormalized continuum 
limit will be considered in the Bethe ansatz framework). Thus we must let — > oo in such a 
way that mo ~ 4a _1 sui7e -2e stays finite. 

The observation of the previous section concerning the doublers provides the basis for the 
following representation of the Fermi fields in the continuum limit 



V>2j ^ Va(xL(ja) + (-YvRija)) 
4>2j+i ^ \fa [XR(ja) + (-) j *7l0'°)) ( 5A ) 



where x an d rj are quantum relativistic Dirac fields. The hopping operator and the local 

charge operator Qj, e.g. for i even, read 

K iii+2 ^ 2a :x[xl- - '• VrVR '■ 

'■Qzj'- - a i : XLXL- + - VrVR- ~{~y (XrVl + VlXr)] ■ 

The symbol : ... : on the r.h.s. refers to the usual normal-ordering for continuum fields in 
the interaction picture. This holds because the operators on the l.h.s. have vanishing vacuum 
expectation value. The complementary cases, namely Q2J-1 and K^j-x^j+l > can be handled 
analogously simply by exchanging right and left modes. These are all the calculations needed 
to obtain the continuum limit of H4, since all quartic terms except the first and the last are 
suppressed as a — > 0. As for the quadratic piece H2, the typical calculation reads 



V4j^2j+2 ^ a Xl( x )xl{x + a) - if R (x)Tj R (x + a) 



+ ( 



Vr( x )Xl(x + a) - x l {x)tjr{x + 



so that, dropping the oscillating terms and developping to first order in a we can calculate the 
non- vanishing terms in the quadratic Hamiltonian as 



V4-V>2j+2 - ^lj + %i>2j - 2a{x ] L d x XL ~ rj R d x r]R) 
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Thus, taking into account the normal-ordering and dropping as above all fast oscillating terms, 
the continuum form of the Hamiltonian reads (here G — > oo as a — > so that the bare mass mo 
stays fixed) 

H — > Hq + H m + Hi n t (5.2) 

where Hq is the kinetic energy 

#o = -i J dx (xrQxXR ~ \\dxXL + rjR d xVR ~ Vl^Vl) , (5-3) 
H m is the mass term 

H m = m J dx (xlXR + XrXl) (5.4) 
and -Hint the quartic interaction 

H int = 2g> J dx (Jpl - Jpl - Jpl - Jpl) . (5.5) 
Here g' = —2 cot 7 and the J's are free-field normal-ordered U(l) currents: 

J%=--xlxa--, J%=-vlva--, a = R,L . 

The nice feature of this result is that all terms surviving the naive continuum limit are manifestly 
Lorentz-invariant, unlike those obtained in the analogous treatment of the XXZ spin chain in 
p). It is natural to regard the mass term H m as 'bare' perturbation of the CFT defined by 
Ho+Hi n t. The troublesome aspect is that 77, the field describing the doublers, does not decouples 
from the putative massive Thirring field x an d prevents a straightforward identification of the 
CFT. 

In order to find the right decoupled description, we use abelian bosonization: 

Xa = M 1/2 :exp(iay/4nu a ): 

1/2 t- nr \ ( 5 - 6 ) 
Va = fi ' :exp(iaV47r?j a ): 

where a = rb (+ = R and — = L), [i is a normalization mass scale, and the fields 

u a (x)= dyJ*(y), v n a {x) = / dy J%{y) 



can be identified with the chiral components of two free massless Bose fields. The symbols : . . . : 
in eq.( |5.6| ) now stand for bosonic free-field normal ordering at the mass scale /j, so that the 
expressions (|5.6|) are effectively /j-independent. 

With the standard rules of abelian bosonization, the Hamiltonian now takes the form, up 
to irrelevant constants, 

h = f dx {{d x u R f + (d xUL f + (d xVR f + (d xVL ) 2 

+ g' [(d x u R )(d x u L ) - (d x v R )(d x v L ) - (d x u R )(d x v L ) - {d x v R ){d x u L )] 

+ m 0/ u :cos[v / 47t(uh + u L )] :| (5.7) 
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Notice that only one boson field is involved in the sine-Gordon interaction, but the mixed 
terms in the third line still couple the two boson fields. In order to elimate them we can use the 
canonical transformations that leave invariant the commmutation rules between left and right 
components of the boson fields: 

[u a (x),u a '(x')] = la6 aa >e(x- x') 
[v a (x),v a >(x')] = \a8 aal e{x - x') 

[u a (x),v a ,(x')] = , (5.8) 

Thus it must be a 0(2,2) trasformation. We find it combining two canonical [7(1,1) transfor- 
mation and a canonical ortogonal SO (2) x SO (2) transformation acting on right and left sectors 
separately. We obtain in this way: 



(5.9) 



(4>r\ 


/ 


r t 


—s 


-t\ 


( UR \ 


Hr 




t —r 


t 


—s 


VR 






-s -t 


r 


t 


UL 




V 


t -s 


t 


-r) 


\vl) 



with 



cosh v cosh A + sinh v sinh A 



s = 



V 7 cosh 2v 
cosh 2 v sinh A + sinh 2 v cosh A 

\J cosh 2v 
cosh v sinh v (cosh A — sinh A) 



v 7 cosh 2v 



and 



tanh(2A) 



2vr 



tanh(2zy) = - sinh(2A) 



In terms of the new fields the Hamiltonian reads 



H 



2vr 2 



1/2 



dx {(d x 4> R ) 2 + (A0l) 2 + m R f + (d^L?) 
+ m fjL J dx :cos[v / 47re A (0/j + 4>l)\ ■ 



(5.10) 



and we see that it correspond to a sine-Gordon model plus a decoupled free massless field. More 
precisely, in passing to lagrangian form, we should scale the fields <j> and £ so that the kinetic 
term is properly normalized. In this way one arrives at the Lagrangian: 



C = \{d^4>) 2 + l{d^) 2 + m 0/ ucos, 



(5.11) 



where 4> = 4>L + 4>R an d £ = £l + £r. The relation of Coleman's coupling constant (5 with g' 
reads 

P 2 1 + ^/tt 



4vr 1 - 2( 5 '/7r) 2 ' 



(5.12) 
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We could now perform the inverse bosonization trick on (ft and £, according to the standard 
rules 0], or with canonical transformation analogous to those done above. This yields at the 
end two decoupled Thirring models, one massive, with Dirac field ift, and one massless, with 
field ip': 

C = tft(i^8^ - m )V> + ^(VW) 2 + W + ^W) 2 • (5-13) 

We see therefore the fermion doubling, charateristic of any local lattice regularization with local 
chiral currents, is completely harmless in our case: it only adds a decoupled massless field to 
the Lagrangian. 

Our derivation is now complete: we have shown that the fermion Hamiltonian (|4.3| ) provides 
a local lattice regularization of the massive Thirring model. The important point is that this 
Hamiltonian is completely integrable, being just the first of an infinite hierarchy of conserved 
charges in involution. One may regard all terms in the lattice Hamiltonian which are of order 
a as irrelevant operators needed to preserve the integrability on the lattice. 

Of course we have performed a 'bare' continuum limit which does not take into account 
renormalization effects. However, the integrability of the model allows to include them exactly 
through the explicit diagonalization of the lattice Hamiltonian. This is carried through by 
means of the algebraic Bethe ansatz, or Quantum Inverse Scattering Method, whose main steps 
will be outlined in the next section. 

On the basis of such Bethe ansatz, on-shell solution, the simultaneous presence of massive 
and massless particles in the relativistic QFT describing the continuum limit was already put 
forward in ref.Q. Here we have related this fact to the well-known phenomenon of fermion 
doubling and have in particular clarified its non trivial off-shell extension. 



6 Main results of the Bethe ansatz 



By definition, the algebraic Bethe ansatz will work in the fermionic formulation just like in the 
standard spin framework. All changes of sign due to the fermionic commutation rules ( |3.9| ) can 
be easily traced down. Wee need not repeat here any derivation, referring to the various review 
articles on the subject (see for instance |J). 

The eigenvectors of the alternating tranfer matrix are written (see eq.( 



\q>) = B(X 1 + i 7 /2| 0)5(A 2 + i-y/2\ 0) . . . B(X r + » 7 /2| 0) 



(6.1) 



where \Q) is the bare vacuum state and the parameters Ai, A2, ■ ■ • , A r satisfy the Bethe ansatz 
equations (BAE) 



sinh(A m + + 17/2) 



sinh(A m + — ij/2) 



N 



sinh(A m - + i-y/2) 
sinh(A m — — ry/2) 



N 



(-) r+1 TT sinh(A m - A n + ij) ^ 
^Jj sinh(A m - \ n -il) ' 
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The eigenvalues read 



A = A A + A D 



a A _ tt sinh(i 7 /2 + A — X k ) 

~ 11 sinh (, 7 / 2 - A + Afc) [b - 6) 

A D = { _y b{X + Q) N b[x _ Q) N TT sinh(3i 7 /2-A + A fc ) 

^ smh(-z 7 /2 + A - A fc ) 



Since one easily verifies that [Q , B] = —B, the BA states ( |6.1[) contain exactly r bare particles. 
Notice also that eigenvectors and eigenvalues depend on 6 both explicitly and through the 
dependence forced on the numbers A& by the BAE. We do not need to consider states with 
more than N bare particles, since they are obtained by particle-hole symmetry, i.e. ipj ipj, 



corresponding to spin inversion in spin language, from the states ( p.l 
As usual, we introduce the so-called counting function [0] 



r 



Z N (X) = N [<f>(\ + 6, 7 /2) + 0(A - 6, 7 /2)] - 0(A - A fc , 7) 

sinh(ix + A) 



k=l 

where 

4>(X, x) = i lo, 



sinh(ix — A) 

has the cut structure chosen so that it is analytic in the strip |9?A| < x. The BAE now read 

Zn(^) = 2n I j , j = 1,2, ... ,r 

where the quantum numbers Ij are always half-odd-integers (we choose to be even). This 
should be compared with the spin formulation where the Ij are half-odd-integers for even r 
and integers for odd r. This appears very natural if we compare the b.c. of antiperiodic fermion 
fields, eq.( 3.12j ), with those corresponding to periodic spins, eq.(pV7|). This difference will play 



a crucial role in determining the U(l) charge of the physical particles. 

The energy (both local and nonlocal) and momentum of a given BA state are calculated 
from eqs. ( grp , ( ggjD , flQjp and (|J). The momentum reads 



P = a' 1 J2 + Aj, 7 /2) - 4>(@ - A,-, 7/2)] (6.5) 



while the local energy is 

. . J_ r h h i 

(6.6) 



A^(e + Ai>7 /2) + A^rt_ A/ .-/2) 



The physical vacuum state, or filled Dirac sea, is the ground state of the local Hamiltonian H, 
that is the lowest possible value of E for fixed N. It corresponds to the unique solution of the 
BAE with N real roots. In the limit — ► oo at fixed lattice spacing a (hence in the infinite 
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volume limit), this solution is described by a smooth density ||. The same applies to all particle 
states characterized by a finite number of holes in the ground state distribution. The energy 
and momentum of one of this holes (a physical fermion) can be calculated exactly to be 

771/ s i f s a . /sinh7rvV7\ (R _v 

E ^ = ~^ 1^-' pM = 2-arctan ^ osh7re/7 j (6.7) 

where tp is the position of the hole in the Dirac sea. This is all rather standard. The important 
novelty concerns the U(l) charge of the holes. 

In the usual spin formulation with periodic b.c, to the removal of a single BA root there 
corresponds the appearence of two holes. Therefore each hole has a renormalized charge Q = 
—1/2. This is clearly incompatible with the interpretation of such holes as fermions, since they 
would not be interpolated by the fermi fields ip n . The sign differences proper of the fermionic 
framework, and in particular the factor (— l) r in eq.( |6.2[) , exactly remedy this. An accurate 
analysis of the phase space available for N — 1 BA roots, using the asymptotic value of the 
counting function Z^{X), shows that only one hole is present in the Dirac sea. This is the 
dressed antiparticle of the original fermion and has charge Q = — 1. As a matter of fact one can 
consider also states with N + 1 BA roots, one of which has imaginary part equal to iir/2: one 
finds the same energy-momentum spectrum of eq.( |6.7[ ), while evidently Q = 1. The dressed 
particle is obtained by particle-hole symmetry. 

The states with one particle and one hole are obtained by removing one real BA root and 
introducing a root with immaginary part equal to 7r/2.This naturally follows by looking at the 
dependence of energy and momentum on the "lattice rapidities" Xj in eqs.(|6.5|) and ( |6.6| ): the 
replacement Xj — * Xj + in/ 2 exchanges the two energy branches in eq.( [4.5[) . 

It is possible to identify the solutions of the BAE corresponding to states with arbitrary 
many fermions and antifermions as well as with breathers (fermion-antifermion bound states 
in the attractive regime 7 > vr/2). A complete and detailed analysis is still lacking in the 
literature (parts can be found in the early BA approaches to the continuum massive Thirring 
model (lq|pq| and in the general study of the BA equations for the XXZ chain jL7|), but is 
outside the scopes of this work. 

For our purposes, it is enough here to examine the continuum limit of the massive part of 

the renormalized energy-momentum. As a , at — ► and — > 00 we find from eqs.( |6.7D the 

Trip 

relativistic expressions in terms of the rapidity 9 = — : 

7 

E = mc 2 cosh 9 , p = mc sinh 9 (6-8) 

provided we identify the renormalized velocity of light and mass as: 

c= — , mc= -e~**h . (6.9) 
27a* a t 

Notice that the velocity of light undergoes a finite renormalization from the bare value cq = a /at 
found before. Of course, we could set the conventional c = 1 by adjusting a/at to 2j/ir. Notice 
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also that eliminating <p from the two relations in ( |6.7D one obtains, still on the lattice 

E = aT x c tanh (^j sin 2 pa + \(mca t ) 2 {l + cospa) 1/2 (6.10) 

with the more precise mass definition meat = 2/ sinh(-7r0/7). This renormalized dispersion 
relation should be compared with the perturbative one, eq.(4.7): apart from an overall factor 
which tends to 1 as G — > oo, all renormalization effects are concentrated in the rescalings 
niQ — > m and Co — > c. In particular the exact spectrum ( |6.10| ) has the same fermion doublers of 
the perturbative one: such doublers are still coupled to the massive modes, as could be checked 
with the direct calculation of the relevant scattering phase shifts. In the continuum limits the 
characteristic momenta of the massless and massive modes get separated by a quantity of order 
a -1 and these scattering phase shifts tend to non-trivial constants. The off-shell decoupling 
shown in the previous section ensures that a proper additional dressing of the massive and 
massless particles exists that decouples them altogether. 



7 On the relation between the coupling costants 



We may now investigate more in details the connection between the parameters of the lat- 
tice Hamiltonian and those of the continuum ones, either bosonic (sine-Gordon) or fermionic 
(massive Thirring). 

The standard relation between (3 in eq.( 5.11| ) and the coupling constant g in eq.( |5.13| ) reads 
(see e.g. Q[|) 

(3 2 _ 1 
Air 1 — g/n 

and differs from the relation ( 5.12j ) derived above with g' . The sine-Gordon coupling constant 
(3 is a regularization-independent parameter, since, for j3 2 < 8tt 2 , the sine-Gordon model can 



be uniquely defined as a perturbated conformal theory [12]. Hence we may safely take /3 as a 
reference parameter to relate g and g': 



9 = 9 



j l + 2g'/n 
1 + g'/ir 



g' i + 



1 + g'/ir 



n-1 



n=l 



(7.2) 



They differ by a formal power series redefinition, as to be expected in the Thirring model, since 
the current -current coupling in two dimensions is cutoff independent but regularization-scheme 
dependent. 



Let us observe, moreover, that the relation (7.2) holds in the interaction picture, since we 
are applying to the interacting sine-Gordon field theory the bosonization rules proper of the 
free bose field. We can relate more precisely (3 to the well-defined lattice parameter 7, and 
then to g' = —2 cot 7 (see eq.(fT5|)), by using exact scaling arguments as follows. The scaling 
dimension of cos (3(f) is P 2 /4tt, since it is fixed by the ultraviolet fixed point, namely the free 
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massless bose field. Through bosonisation cos (3(f) maps into ipip, which enters the lagrangian 
of the massive Thirring model multiplied by mo- Hence mo has scaling dimension 2 — 1 /A-k. 
From the exact Bethe ansatz solution one learns that the physical mass scale is proportional to 



exp(— 7T0/7) (see eq.(|6.9P). On the other hand, eq.(4.8) shows that mo scales like exp(— 20), 
so that it has scale dimension 27/77. Therefore we must have 27/71 = 2 — /3 2 /4-7r, which is the 
exact relation we sought. 

This argument is quick but rather too sketchy. A more precise derivation goes at follows. The 
redefinitions of the normalization mass scale \x and those of the the bare mass mo are connected 
by the normal-ordering renormalization group [11 1, in order to keep mvipip = mo/J : cos f3(j) : 
invariant. This leads to the relation 

mon ( fi x A 



m /x' 



where A = /3 2 /47r. On dimensional grounds, the physical mass scale has the form 

m = m f{p,z) , z = — 



/i 



and must be renormalization-group invariant, that is 



m = m A A +7(/^A A + 2 ) , X = ±. 

M 

Hence / is a homogeneous function of z and we obtain 

m = m z- y /(/3,l)=m^V/(/3,l) (7.3) 

where y = A -±±. 

The exact Bethe ansatz solution of the lattice model provides the following relation for the 
fermion mass in the continuum limit 

m a A*U^>h = ^( a -if (^lY' 21 (^S" 1 " 1 

a 2 7T 7r V a / V sin 7/ \ 4 / 

where eq. (|4.8|) was used in the second equality. Choosing /x = a^ 1 and a /at = 27/71, to enforce 
c = 1, we obtain, comparing to eq.(7.3), 



4vr/(/J 2 -87r) 



The relation between 7 and (3 is that we found above. In addition we found an expression for 
/(/?, 1). Of course this expression is scheme-dependent. 

Coming back to the Thirring coupling constants g and g', we have the following situation: g is 
defined through bosonization of the massive Thirring model alone and is given by g = ir— 4-7T 2 //? 2 
(see eq.( [7.1[) ). Hence we have the exact relation 

(7-4) 

2 7T — 7 
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g' may be analogously defined through bosonization of the complete continuum hamiltonian 
(|5.2j ) (which contains the fermion doublers). This leads to the exact relation (|7.2| ). On the 
other hand we have the relation g' = —2 cot 7 (see below eq.( |5.5[) ), which follows from the 
continuum limit of the lattice hamiltonian in the interaction picture (using free-field normal- 
ordering), and therefore is only approximate or "bare". Combining eqs.( [7.2[ ) and Q7.4| ) we obtain 
the exact relation 

2g'l + 2g'/ir _ 1-27/vr 
7r 1 + g'/ir 1 — 7/7T 

which can be regarded as the renormalization of the bare relation g' = —2 cot 7. 



8 Final comments and outlook 

The local lattice regularization of the massive Thirring model presented here applies equally 
well to the vast class of integrable models already under control by means of the standard light- 
cone approach. The local character of the lattice Hamiltonian should help in the field-theoretic 
understanding of these models, since it allows for a better control of the continuum limit and 
an easier identification of each model as a perturbed CFT. From this field-theoretic point of 
view, the most important step remains the proper definition of the local lattice fields in terms of 
which the 72-matrices are to be written. When this is done, the Hamiltonian as well as all other 
conserved charge, either local or nonlocal, would follow by the standard techniques of vertex 
models, since only the algebraic properties of the .R-matrices and the permutation operators 
are needed. In the case of the massive Thirring models this program hass been here pursued 
explicitly starting from the local .R-matrices written in terms of canonical lattice fermi fields 
(eq.^3) and handling the continuum limit as in section |H[ 
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